Experimental realization of SQUIDs with topological insulator junctions 
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Abstract 

We demonstrate topological insulator (Bi2Te3) dc SQUIDs, based on superconducting Nb leads coupled 
to nano-fabricated Nb-Bi2Te3-Nb Josephson junctions. The high reproducibility and controllability of the 
fabrication process allows the creation of dc SQUIDs with parameters that are in agreement with design 
values. Clear critical current modulation of both the junctions and the SQUID with applied magnetic fields 
have been observed. We show that the SQUIDs have a periodicity in the voltage-flux characteristic of <I>o, 
of relevance to the ongoing pursuit of realizing interferometers for the detection of Majorana fermions in 
superconductor- topological insulator structures. 
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A three-dimensional topological insulator is an insulator in the bulk, but has conducting sur- 
face states that can be described by means of a Dirac cone in which the spin is locked to the elec- 
tron momentum [iJlO]. Electrons with opposite spin have opposite momenta, which suppresses 
backscattering, rendering these materials interesting for low power electronic devices. Combining 
the helical Dirac fermions with a superconductor [111- 13] may lead to the artificial creation of the 



elusive Majorana fermion lll4ll . In the search for the Majorana fermion, efforts have been made 



to contact a topologica 



have been reported nl5 



insulator (TI) to a superconductor (S). Supercurrents in S-TI-S junctions 
- llTII . The combined evidence for a ballistic Josephson supercurrent and 



the presence of topological surface states, show that topological Josephson junctions can be made, 
despite the presence of a conductivity shunt through the bulk of the TI crystal [[itI] . 

In this Letter, we demonstrate dc SQUIDs consisting of Josephson junctions with topological 
insulator surface states as barrier layer. These interference devices potentially serve as a basis to 
detect the Majorana fermion. Majorana bound states may appear in the vortex of a topological 



superconductor and at S-TI interfaces 



131] . In the latter, time reversal symmetry in the topo- 



logical insulator needs to be broken, for example by incorporating a magnetic insulator layer or by 
applying an external magnetic field. Intriguing devices are proposed to identify the appearance of 



this exotic particle 111 iM 131.1 1811. Since the Majorana fermion is charge neutral and is a zero energy 



state, most proposals rely upon quantum interference devices. For example, a current-phase rela- 
tionship with a An periodicity 111 ill might result from the interplay between the Majorana fermion 
and a superconductor. We therefore study the current-phase relationship of the topological dc 
SQUID. 

Figure \T\ shows two superconductor Nb - topological insulator Bi2Te3 dc SQUIDs fabricated 
on one Bi2Te3 flake. We have fabricated poly crystalline Bi2Te3 samples with a common c-axis 
orientation using the Czochralski method as described elsewhere lll9ll . Using mechanical exfo- 
liation, Bi2Te3 flakes ranging from 30 nm to 1 |jm are transferred to a Si substrate. The Bi2Te3 
flakes are smooth on the nm scale over areas of several \im^ scale. Figure [2] shows an Atomic 
Force Microscopy image of the surface of a typical Bi2Te3 flake, revealing the 1 .0 nm quintuple 
unit cell layers of Bi2Te3. After exfoliation, a superconducting (Tc = 9K) Nb-layer (200nm) is 
sputter deposited with a 5nm Pd layer deposited in situ on top to protect the Nb against oxidation. 
Electrodes are defined by photolithography. The Nb on top of the Bi2Te3 flake makes a strong su- 
perconducting contact with the Nb on the substrate. Finally, nanojunctions are defined by lift-off 
e-beam lithography and sputter deposition of Nb. The substrate is slightly conducting (p = 5ficm) 
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FIG. 1: Scanning Electron Microscopy image of topological insulator nano dc SQUIDs. Two SQUIDs with 
length ^=120nm are designed on an exfoliated Bi2Te3 (200nm height) flake. The two Nb arms are differently 
sized resulting in an asymmetric SQUID. Nb is sputter deposited and defined by e-beam lithography. Inset 
shows a Josephson junction. 



FIG. 2: Atomic force microscopy image of an exfoliated Bi2Te3 surface. The step edges are 1.0 nm high, 
corresponding to the Bi2Te3 quintuple unit cell. These nanometer flat surfaces span an area up to 50 x 
50 laml 

at room temperature to increase the resolution of e-beam lithography, but is completely insulating 
at low temperatures. Prior to deposition, in situ Ar-ion etching is performed in order to make trans- 
parent contacts. It is found that Ar etching roughens the Bi2Te3 surface, probably by preferential 
etching, but 2 minute etching at 50 eV leaves a Bi2Te3 surface with 1-2 nm roughness and a trans- 
parent contact. Throughout the fabrication process the Bi2Te3 surface in between the Nb leads 
has only been covered by resist, leaving the surface of the junction barrier layer unaffected from 
Ar-ion etching and deposition steps. Two separate electrodes connected by Nb over the Bi2Te3 
flake, for example the positive current and voltage leads, have a superconducting contact with a 
critical current exceeding 30 mA, ensuring a large supercurrent through all layers. 

We have designed two SQUIDs on one Bi2Te3 flake (200nm height), see Fig. \\\ The nanojunc- 
tions are 2 |am wide and have a length / = 120 nrti. The inductance ratio between the two arms is 
about A = 0.2. In a square washer approximation [20] the total effective area of the SQUID is 920 
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FIG. 3: IV -curve of the topological dc SQUID with atypical critical current of 30|aA. The junctions are 
2^m wide and have an electrode separation of 120 nm. The intrinsic bulk shunt has reduced the IcRn 
product to 10 ^V. Inset shows the IF-curve under magnetic field modulation. 

|am^, and the estimated inductance 46 pH. Both SQUIDs showed similar behavior, but in the rest 
of the paper we focus on the left SQUID of Fig. 1. 

Below 6 K, the superconducting proximity effect induces a Josephson supercurrent through the 
junctions and at 1.4 K the dc SQUID has a critical current of 30 |jA, as shown in figure [3l This 
results in a 2D critical current per width of the junction of 7.5 A/m, which is within 10% of the 
individually measured critical current density of junctions on different flakes. The SQUID had 
a critical current constant within 10% over three cooldowns running over several weeks. Due to 
shunting of the junctions results we find a relative low IcRn = 10 |aV. Ballistic junctions with 
^0 ~ 80nm, Tc = 6K and /=120 nm, have an estimated IcRn ~ 130 — 260|aV, but the shunt due 
to the bulk conductance reduces the characteristic voltage to a few percent of the expected value 
Inll . consistent with the observed 10 |aV. In order to obtain higher IcRn values, electrical gating or 
chemical substitution, e.g. Bi2Se3Te3 [|2ll] . could be used. The current-voltage characteristics are 
intrinsically non-hysteretic because of the low capacitance of the lateral geometry of the junctions 
and the high transparency of the interfaces. 

Applying an external magnetic field causes oscillations of the superconducting critical currents 
due to interference of the two arms, shown in|4l This modulation of Ic unequivocally demonstrates 
the correct operation of the Nb-Bi2Te3 SQUID. The voltage modulation is strongest at a bias of 
approximately 40 |aA, which corresponds to a current close to the critical current of the SQUID. 
The critical current modulation is ~ 35% of the total critical current. Using a simple model 
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FIG. 4: Voltage modulation the dc SQUID, (a) The critical current of the SQUID can be fitted (thin line) by 
a model that accounts for the asymmetry of the design, (b) The critical current of the individual junctions is 
suppressed at 350^T (Fraunhofer pattern). SQUID oscillations are still observable. Inset shows the SQUID 
oscillations at small fields (left) and large fields (right) (c) The V-0 relationship reveals an asymmetric 
response as well as a shift for increasing bias currents. The solid lines represent equicurrents for it 30-70 
MA. 

based on a sinusoidal current-phase relation of the junctions, and taking inductance and possible 
asymmetries into account, we can fit the field dependence with high accuracy. The asymmetrical 
inductance causes an asymmetry in the SQUID current-phase relationship, see Fig. HI We also 
observe a shift in the current-phase relationship when increasing the bias current. The two arms 
of the SQUID have different inductances, which results in a growing field threading the SQUID 
for increasing current. The 0.75 $0 shift in the V-0 relation at 90|aA bias can be accounted for by 
a lOpH inductance difference between the two arms. This corresponds to the fitted values for the 
total inductance of the SQUID, = ^ = 55pH and the asymmetry factor A = 0.2, resulting in 
an inductance difference of 1 IpH. At the optimum bias, the sensitivity of the SQUID is 15 |aV/<l>o- 
In the limit /9l(= 5) ^ 1, the sensitivity of a SQUID can be estimated by R/L [|22|l which yields 
12 nV/$o. 

As expected, the SQUID modulation frequency is much larger than that of the individual junc- 
tions, see Fig. HI A SQUID oscillation corresponds to 1.9 |aT, while the critical current of the 
junctions is suppressed at 350 |aT. This corresponds to about 180 SQUID oscillations in a junc- 
tion oscillation, which is slightly lower than expected by comparing the enclosed areas. However, 
for the enclosed area of a junction, the Josephson penetration depth and flux focussing has to be 
included, which increases the junction effective area uM. 



Since the supercurrent is carried by topological surface states [|l7ll we can study the current- 



5 



phase relationship of toplogical junctions by means of our dc SQUID. In a dc SQUID with topo- 
logical insulator surface states as interlayers, the phase difference over the junctions due to a given 
field is equal to a standard SQUID. However, the presence of Majorana fermions can cause a 
Jc = /osin(f/)/2) junction current-phase relation, resulting in a An periodic dependence of the 
junctions nlu\ . To test the periodicity of the SQUID, the magnetic field from the coil is calibrated 
with a Hall sensor and the effective area of the junction, 920|am^, is estimated using the square 
washer approximation ll20ll . A 2% periodic dependence within 7 % accuracy is calculated, enough 
to exclude An periodicity. The absence of a An periodic dependence might result from the absence 
of a ferromagnetic insulator (breaking time reversal symmetry), or due to relaxation to equilibrium 



states as was predicted theoretically 1123 . 



2M. 



In conclusion, we have demonstrated the Nb-Bi2Te3 SQUIDs. As shown previously, the super- 
current is carried by the topological surface states, thereby allowing the study of the current-phase 
relationship of superconductor- topological insulator structures. From the V-0 characteristics we 
deduced a 2n current-phase relationship. The high reproducibility allows to study the noise prop- 
erties of these devices in future experiments and to eventually include magnetic insulators in order 
to break time reversal symmetry to create Majorana bound states. 
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